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We propose fractional Fokker-Planck equation for the kinetic description of 

^D . relaxation and superdiffusion processes in constant magnetic and random elec- 

tric fields. We assume that the random electric field acting on a test charged 
particle is isotropic and possesses non-Gaussian Levy stable statistics. These 
assumptions provide us with a straightforward possibility to consider formation 
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c/2 \ of anomalous stationary states and superdiffusion processes, both properties are 

. ^h ' inherent to strongly non-equilibrium plasmas of solar systems and thermonuclear 

^- devices. We solve fractional kinetic equations, study the properties of the solu- 

^G . tion, and compare analytical results with those of numerical simulation based 

on the solution of the Langevin equations with the noise source having Levy 
stable probability density. We found, in particular, that the stationary states 
K> , are essentially non-Maxwellian ones and, at the diffusion stage of relaxation, the 

characteristic displacement of a particle grows super diffusively with time and is 
inversely proportional to the magnetic field. 
PACS: 05.10 Gg, 05.40. Fb 



1 Introduction 

Anomalous random motions and related transport phenomena are ubiquitous in na- 
ture. In these phenomena the laws of normal diffusion (ordinary Brownian motion) 
are altered, e.g., the mean square no longer increases linearly in time, but instead 
grows slower (sub-diffusion) or faster (super-diffusion) than the linear function. There 
are a lot of examples from very different fields of applications, see reviews (T|], and 
references therein. The anomalous random motions often exhibit long-time and/or 
space dependence as well as multi-scaling, or multifractal, behavior @], 0. These cir- 
cumstances require to go beyond the theory of (relatively simple) random Markovian 
processes as well as beyond the theory of (mono) fractal, or self-affine, processes. The 
systems, in which anomalous random motions occur are usually essentially non-linear 
and, in this sense, the random motions are non-linear ones; this circumstance again 
greatly complicates the problem of an adequate statistical description. 

The two basic anomalous fractal random motions are of particular importance, 
namely, fractional Brownian motion || , and the Levy motion, whose theory has begun 
from the works of P. Levy || . The former motion is characterized by long-range time 
correlations, whereas the latter one is characterized by non-Gaussian statistics ; in 
this case the increments of the process may be independent (Levy stable processes or 
ordinary Levy motion 0) or have an infinite span of interdependence (fractional Levy 
motion) §, [1, . 

The theory of Levy stable distributions and stable processes naturally serves as 
the basis for probabilistic description of the Levy motion, since stable distributions 



obey Generalized Central Limit Theorem, thus generalizing Gaussian distribution [I I 
However, the application of the theory of stable processes is limited because of the 
infiniteness of the mean square and discontinuity of the sample paths. The finite 
sample size and boundary effects play a decisive role, thus modifying stable probability 



laws ("truncated Levy distributions") |12|], and violating the property of self-affinity 
("spurious multi-affinity" ) [|T3| . 

The peculiarity of anomalous random motions is that they cannot be described 
by the standard Fokker-Planck equation because the basic assumptions, namely, the 
Markov property and the local homogeneity of the space do not hold in these cases. 
The use of differential equations with partial fractional derivatives is a perspective way 
for describing such processes. One of the implementations of such an approach is the 
use of different forms of fractional kinetic Fokker-Planck equation or the corresponding 
forms of the Master equation with fractional derivatives. Recently, kinetic equations 
with fractional partial derivatives have attracted attention as a possible tool for the 
description of diffusion and relaxation phenomena, see review and references therein. 
However, fractional calculus is far from being a regular tool in physics community, 
and the solutions to fractional kinetic equations are known in a very few cases. The 
development of the theory requires, from one hand, the development of microscopic 
foundations of fractional kinetics and, from the other hand, the development of powerful 
regular methods for solutions to fractional equations. 

The various processes in space and thermonuclear plasmas could serve as impor- 
tant applications of fractional kinetics. Indeed, many of the current challenges in solar 
system plasmas as well as in plasmas of thermonuclear devices arise from the funda- 
mentally multiscale and nonlinear nature of plasma fluctuation and wave processes. 
Anomalous diffusion and plasma heating, particle acceleration and macroscopic trans- 



fer processes require to go beyond the "traditional" plasma kinetic theory. Fractional 
kinetics can be useful for describing such processes, just as it occurs in other fields 
of applications. Our paper is a step in this direction. We consider the motion of a 
charged Levy particle in a constant external magnetic field and random electric field 
obeying non-Gaussian Levy statistics. Our problem is a natural generalization of the 
classical example of the motion of a charged Brownian particle ||14|| . We solve the frac- 
tional Fokker-Planck equation with fractional velocity derivative, study the relaxation 
processes in phase and real spaces as well, and estimate fractional moments of energy 
and coordinate. We also perform numerical modelling based on the numerical solution 
to the Langevin equations and demonstrate qualitative agreement between analytical 
and numerical results. 

2 Fractional Fokker-Planck equation for charged par- 
ticle in magnetic field 

The history of fractional Fokker-Planck equation (FFPE) for the probability density 
function (PDF) f(f,v,t) in the phase space goes back to the papers by West and 
Seshadri flT5| , and by Peseckis JTBJ. Here we recall briefly the arguments used when 
deriving FFPE. It is well known that usually the derivation of classical kinetic equations 
for the Brownian motion is based on the assumption of the finiteness of the second 
moments of the PDF. This way is not useful here, because, as we shall see, the second 
moments diverge. Thus, it is expedient to explore the method used by Chandrasekhar 
17] for the derivation of the Fokker-Planck equation for the Brownian motion. His 



method does not require the finiteness of the second moments. In fact, for fractional 
kinetics, the modification of Chandrasekhar's method was proposed for the first time in 
Ref. [I6|. Here we proceed mainly along the derivation from Ref. [fTH |. We consider a 



test charged particle with the mass m and the charge e, embedded in constant external 
magnetic field B and subjected to stochastic electric field 8 (t) .We also assume, as 
in the classical problem for the charged Brownian particle fL4|] , that the particle is 
influenced by the linear friction force —vmv, v is the friction coefficient. For this 
particle the Langevin equations of motion are 

dr 

* = v • (21) 

dv e r ^ -, e =♦ 

— = — [v x B \ — vv H 8 . 

at mc ' m 

The statistical properties of the field 8 (t) are assumed to be as follows. 

1. 8 (t)is homogenous and isotropic. 

2. 8 (t) is a stationary white Levy noise. 

The first assumption is the usual one when dealing with the motion of charged 
particle in a random electric field. In subsequent Sections we consider two possibilities: 

(i) 8 (t) is a 2-dimensional (d = 2) isotropic field in the direction perpendicular to 
external magnetic field; in this case the motion along the magnetic field is neglected, 
and 

(ii) 8 (t) is a 3-dimensional (d = 3) field. 

The second assumption has a profound meaning. 

Indeed, if 8 (t) is a white Gaussian noise, then we encounter with a classical Brown- 
ian problem and, by using Eqs. (2.1), we arrive at the Fokker-Planck equation and get, 



as the consequences, Maxwell stationary PDF over the velocity, exponential relaxation 
to the Maxwell PDF and the normal diffusion law for the particle motion in the real 
space, as well. Instead, the non-Gaussian Levy statistics of the random force in the 
Langevin equation (2.1) provides us with a simple and straightforward, at least, from 
the methodical viewpoint, possibility to consider abnormal diffusion and non-Maxwell 
stationary states, both properties are inherent to strongly non-equilibrium plasmas of 
solar system and thermonuclear devices. Returning to Eq. (2.1), it follows from the 
assumption 2 that the process, which is an integral of £ (t) over some time lag At, 

t+At 

L(At) = I dt'£(t') , (2.2) 

t 

is an a-stable isotropic process with stationary independent increments 0, ||, whose 
characteristic function is 

Qj l (k, At) = exp f-D £ k * At) , (2.3) 

where a is called the Levy index, < a < 2 , and the positive parameter Dg has the 
physical meaning of the intensity of the random electric field. If a = 2 , Eq. (2.3) is 
a characteristic function of the Wiener process, and, after applying Chandrasekhar's 
procedure we arrive at the Fokker-Planck equation. But if a < 2 then, by applying 
the procedure described in detail for the one-dimensional case in Ref. |1|, we arrive 
at the fractional Fokker-Planck equation for the charged particle in the magnetic field 
and random electric field: 

% + *% + n & x ^ % = v lfv m ~ D { - A " )a/2 f ' (2 ' 4) 

where £1 = eB/mc , D = e a D£/m a and (— A^) a ' is the fractional Riesz derivative 
over the velocity. This operator is defined through its Fourier transform as 

(- A,) a/2 f (f, v, t) - \k\ a f (x, k, t) (2.5) 

where / is the characteristic function, 

/ f x, fc, t J = ( exp f ixr + ikv J ) , (2.6) 

the brackets (...) imply statistical averaging. 

An explicit representation of the Riesz derivative is realized through hypersingular 



integral, see the monograph |19| containing detailed presentation of the Riesz differen- 
tiation. We also note that at a = 2 Eqs. (2.4), (2.5) are reduced to the Fokker-Planck 
equation of the Brownian motion. In the next Sections we get the solution to Eq. (2.4) 
and consider physical consequences. 

3 Solution to fractional Fokker-Planck equation 

In this Section we solve Eq. (2.4) with initial condition 

f(?,if 1 t = 0) = 6(f- f ) 5 (v - v ) (3.1) 



in the infinite phase space. We pass to the characteristic function (2.6), which obeys 
the equation 



-J- + (-k + n (k x b) + uk) -L = -D 
8t \ \ ) ) dk 
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with the initial condition 



/ f x, k, t — J = exp f ixr + ikv ) . 



(3.2) 



(3.3) 



Equation (3.2) can be solved by the method of characteristics. The equations of char- 
acteristics are 



— = vk + Q[kxb) — h 
at 



(3.4a) 



dt 



/■ 



(3.4b) 



Denote k! (t') the value of k at time instant t'. Then, the solution to Eq. (3.4a) can be 
found after lengthy, but straightforward, calculations: 

K = e^*"*') | (K'b) b + bx (k' xb) cos Q(t- f) 



+ (it' x b)smtt(t-t') 



(3.5a) 



it' = e v{t '- i] { (itb) b + bx (it x b) cos Q (t' - t) 



+ [K x b)smQ(t' -t) 



(3.5b) 



where K = k-G,K' = k'-G , 



G = A_Z + JL [ i, x ( 2 x f» 






b x x 



(3.6) 



Now, from Eq. (3.4b) we can get expression for the solution of Eqs. (3.2), (3.3): 



/ ( x, k,t) = f ( x, k' {t' = 0) , t = ) exp <^ -D / di' fc' (t 



(3.7) 



In Eq. (3.7) k' {t') is expressed through k with the use of Eqs. (3.5), (3.6). In the next 
Sections we consider the peculiarities of the relaxation process and of stationary states 
realized in the framework of this solution. 



4 Homogeneous relaxation and stationary states 

In this Section we consider homogeneous relaxation, d/df = in Eq. (2.4). Obviously, 
it corresponds to the particular case x = in the equations of Section 3. Setting x = 
in Eqs. (3.5)-(3.8) we get 



f(k,t) =f x = 0,M =exp\ik v (l-e-°^) 
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fc&) 6+ lb x fex6 cos fit + (6 x k) sin fit 



(4.1) 



(4.2) 



(4.3) 



Let us consider the case Vq = for the sake of simplicity. The relaxation process, 
whose characteristic function is given by Eq. (4.1), is not an a-stable process with 
independent increments, since D^ is not a linear function of time, see Eq. (4.2). The 
stable process arises asymptotically at small times only, 



t « t v = 1/au, 



(4.4) 



when the exponent in Eq. (4.2) can be expanded into power series. On the other 
hand, after exponential relaxation to the stationary state, that is, at t » t v the 
stochastic process v (t) becomes asymptotically stationary process with the stable PDF 
independent of t, 



fstlk 



D 
exp [ 

av 



(4.5) 



We also note that stationary PDF does not depend on the magnetic field. 

Another interesting point is related to stationary solutions of fractional kinetic 
equations. In the theory of Brownian motion equilibrium Maxwell PDF over velocity is 
reached at t » \jv. It is characterized by the temperature T of surrounding medium. 
The following relation exists between the parameter D and the friction coefficient v : 



D 



vksT 



m 



(4.6) 



where &b is the Boltzmann constant. The temperature T is a measure of a mean kinetic 
energy of the Brownian particle: 



(E) 



m (v 2 ) ksT 



(4.7) 



Equations (4.6) and (4.7) are examples of fluctuation-dissipation relations. For this 
case the source in the Langevin equation is called the source of internal fluctuations. 
Relations (4.6) and (4.7) may not be fulfilled, as it takes place, e.g., in auto-oscillation 
systems [|(J . In such a case one says that there is the source of external (relatively to 
the system considered) fluctuations in Eq. (2.2). However, Maxwell exponential form 
of stationary solutions retains. As to the Levy motion, the fluctuation - dissipation 
relations can not be fulfilled, at least, because of the infinity of square velocity: (v 2 ) = 
oo for < a < 2. Therefore, we can only say about the Langevin source as about 



the source of external fluctuations. Moreover, the stationary solutions do not possess 
Maxwell form but instead more general form of stable distributions. We also note that 
at present there is no theory of equilibrium state basing on stable PDFs. 
We further study energy distribution in the stationary state, 



fst(E) = Jdvf st {^)8\E 



mv 



(4- 



dvSlE 



exp l-ikv) fst(k) . 



mv 2 \ T dk 

~)J (^r 

We recall, that the two possibilities can be considered, namely, (i) the random 
electric field is isotropic in the plane perpendicular to the external magnetic field, and 
(ii) the field is isotropic in the three-dimensional space. In the former case k and v 
are two-dimensional vectors in Eq. (4.8), d = 2, whereas in the latter case k and v are 
three-dimensional vectors, d — 3. We consider both cases. After some transforms we 
get from Eq. (4.8), 



fst (E) 



fst (E) 



— f dk-kjA k\— | exp (-Vk a ) , d = 2, 
m J \ V m I 

o \ / 



nm 



dk ■ k sin k 



2E 



in 



exp(-Vk a ), d = 3, 



(4.9a) 



(4.9b) 



where T> = D/av. The integrals in Eq. (4.9) can be easily calculated in two particular 

cases: 

1. a = 2 



fst (E) 
Lt (E) = - 



1 



2mV 
1 



exp 



2mV 



d = 2, 



¥(mD) 



3/2 



exp 



2mV J ' 



(4.10a) 
(4.10b) 



which are the well-known results of the theory of Brownian motion |20 
2. a = l. 



fst (E) 



1 /mP 2 



2\ 1/2 



fst (E) 



(E + mV 2 /2) 3/2 
V2mV 2 ^E 



^2 ' 



d 



(4.11a) 



(4.11b) 



vr (E + mV 2 /2Y 

Since the a-stable distribution with a = 1 is called the Cauchy distribution, Eq. (4.11) 
corresponds to the case of the Cauchy motion. 

From Eq. (4.9) it follows that at large energies f s t (E) has a power law asymptotics 
for < a < 2, 



f st (k) oc E-^V ; 



(4.12) 



and, thus only the moments of the order g less than a/2 are finite for a < 2. For the 
moments of the energy, 



{E") = 1 dEE"f st (E) (4.13a) 

o 



we get 



(&) = (2mf D 2 '/«r (1 + q) F [V 2q/ ^ , d = 2, 

r(i-g) 



(4.13b) 



(E*) = A (2myV^ . **** r (g ^j 3 { 2 + g) , d = 3, (4.13c) 

x 7 v 7 ^ sm2vrg/a T(2q/a) K J 

where g < a/2 < 1. The particular cases a = 2 and a = 1 can be also obtained from 
Eqs. (4.13) or by direct using f st (E) from Eqs. (4.10), (4.11): 
a = 2. 

(E q ) = (2mV) g T (1 + q) , d = 2, (4.14a) 

( ^ H M^ r (, + 3\ ,^ 3 , ,4,4b) 



for all g > 0. 
a — 1. 

< £« >= f^-) -1= r (1 + g) r fl - g 1 , d = 2, (4.15a) 



7T 



< B ">= i zirf r (S + 9]r(i-9), rf = 3, (415b) 



( -^r(? + ,)r(i 

for g < 1/2. 

We carry out numerical simulation based on the solution of the Langevin equations 
(2.1) with a two-dimensional isotropic white Levy noise S(t). The case of a strong 
magnetic field is simulated. In Fig. 1 typical dependencies E (t) (t is a discrete time, 
At = 10~ 3 is a time step) are shown on the left for a) a — 1.95 and c) a = 1.1. 
With Levy index decreasing the "jumps" on the trajectories, or "Levy flights", become 
larger; this effect is due to the power law tails of the PDF of the white Levy noise in 
the Langevin equation. In this Figure and below, in Figs. 2-4, the parameters used 
in simulation are Q = 2, v — 0.07, and D — 1. At the right, Figs, b) and d), the "Levy 
flights" are shown on the (v x ,v y ) plane. Again, large "jumps" are clearly seen in the 
bottom figure. 

In Fig. 2 the stationary PDFs f st (E) are shown for a two-dimensional problem at 
the left, Figs, a), c) and for a three-dimensional problem at the right, Figs, b), d), 
respectively. At the top, in Fig. a), b), the linear scale is used, whereas in the bottom, 
in Figs, c), d), the PDFs are shown in the log-log scale. 

The PDFs estimated according to Eqs. (4.9) for the Levy index a — 1.1 are depicted 
by solid lines, whereas the PDFs estimated according to Eqs. (4.10) for a = 2.0 
(Brownian motion) are depicted by dotted lines. The bottom figures clearly show the 
power asymptotics of the PDFs. The black points on the left figures indicate the PDF 
obtained in numerical simulation of a two-dimensional problem, the parameters used 
in numerical simulation are the same as in Fig. 1. The quantitative agreement between 
analytical and numerical results is obvious. 

In Fig. 3, as the result of numerical solution of the Langevin equations, the moments 
(E q ) of the energy are shown for the Levy index a = 1.6 and for different orders q, see 



from the bottom to the top: q = 0.12, which is less than a/2, q = 0.8, which is equal 
a/2, and q = 2.0, which is greater than a/2. The stationary level of the q-th moment, 
estimated according to Eq. (4.13b) is indicated by dotted line in the bottom figure. It 
is seen that, with q increasing, the moments of the energy strongly fluctuate; this is 
the numerical manifestation of the fact that these moments diverge at q > a/2. 



5 Non- Homogeneous Relaxation and Superdiffusion 



We turn to the relaxation in non-homogeneous case. Since general analysis of Eqs. 
(3.5)-(3.8) is rather complicated and taking in mind that we already have information 
about velocity relaxation, we study evolution of a simpler PDF instead of f(r,v,t), 
namely, 

/ (f, t \f , vo) = dvf (f, v, t \f , vq ) , (5.1) 

whose characteristic function is 

/ (x,t\f ,v ) = f (x,k = 0,t\r Q ,v o ) , (5.2) 

and the characteristic function in the r.h.s. of Eq. (5.2) is given by Eqs. (3.5)-(3.8). 
Putting k = in Eqs. (3.5), we get 



k'{t' 



Kb\b 
v 



1-e 



u(t'-t) 



+ 



cos Q it' - 1) + — sin tt (£' - t) 
v 



+^ (b x x) { 1 - e^-V [cos n(t'-t)-^- sin Q (t' - £)] } 



(5.3) 



In the absence of the magnetic field, B = , we get, using Eqs. (5.3), (5.2)and 
(3.7), 



f{x,t \ro,v ) = expOxr H (1-e ) — 



vt\ a 



dr (1 - e~ ut ) 



(5.4) 



For one- dimensional case this result was obtained in 18 



For the case of a strong magnetic field, fi >> v , we get, again by using Eqs. (5.3), 
(5.2) and (3.7), 



/ (x, t |r , vq ) = exp ^ ixf + ik v — D dt' k! (£' 

o 



(5.5) 



where 



k'(t') = V J (l-e^ r ) + — ffrx (xxfrjjsinfir + (5.6) 

1 — e~ UT cos fir (? \ . 

+ s (6xx), T=t-f, 

and fc is given by Eq. (5.6) at t' — (that is, t = t) . 

For simplicity we put r = , and fJo — , as in previous Section. Further, we are 

interested in diffusion across the magnetic field and, therefore we set ( Kb 1 =0 in Eq. 

(5.6), thus not considering the motion of a particle along the magnetic field. For the 
characteristic function we get 



/ (x, t) = f (x, 1 10, 0) = exp < —^x a f dr (l - 2e~ VT cos fir + e - 2uT ) a/2 } . (0.7) 

— * 

where x is two-dimensional vector in the plane perpendicular to B , x = |x| . 

At r >> 1/V (the diffusion stage of relaxation) the expression (5.7) gets a simple 
form, 

/(x,t) = exp|-^x«t} , (5.8) 

which is the characteristic function of an a-stable isotropic process, compare with Eq. 
(2.3). Now we consider the PDF and diffusion process in more detail. By taking inverse 
Fourier transform from Eq. (5.8) we get the PDF, 



t 

f (r, t) = j J 0e—7 (*, t)=J ^ Jo (xr) exp (-^*°) , r=\¥\. (5.9) 

o 

The particular cases of Eq. (5.9) are 

1. a = 2, 

/(r '^4^ eXP (-^ r J' (5 - 10) 

and 

2. a = 1, 

/(r,t)= - Dt/ " . /2 • (5.11) 

2vr (r 2 + (/Jt/fi) 2 ) 7 

In general case < a < 2 the asymptotics of the PDF at large r behaves as 

f(r,t)oc-^— . (5.12) 



It implies that the q-ih moments of r diverge at q > a . The expression for the 
moments follows from Eq. (5.9): 



(r q ) = M r (£; q, a) 



(Dt) 



l/a 



Q 



C 2 (q; a) 



(5.13) 



where 



OO CO 

C 2 (q; a) = dx X K+ q I dn 2 x 2 J (xi>c 2 ) e~^ ? ■ (5.14) 



The integral over x 2 behaves as yc v 2 ~ a at large H\ , thus C 2 (q; a) diverges at upper 
limit of K\ at q > a. The particular cases following from Eqs. (5.13), (5.14) are as 
follows: 



1. a = 2, 



K^f V2r(i+ f>' «>» 



(5.15) 



and 



2. a = 1 



1 (Dt\ q m f q 1 q 

N ' 2\n J V 2 



< q< 1 



(5.16) 



where B is the beta-function. 

From Eq.(5.15) we get a classical diffusion law for the square displacement of 
charged particle across the magnetic field: 



<r>> 



cx 



t 



(5.17) 



One can introduce "the typical displacement" of a particle defined as 

Ar=(r q ) 1/q . (5.18) 

From Eq. (5.15) it follows that at any q (not only at q — 2) we have for the Brownian 
particle, 



Ar a 



t l/2 

IT 



(5.19) 



with the prefactor which, of course, depends on q. We recall that usually just t and 
(especially in plasma physics) -B-dependences serve as indicator of normal or anoma- 
lous diffusion. Therefore, we can use the typical displacement (5.18) as a measure of 
anomalous diffusion rate at < a < 2 and any q < a. Indeed, it follows from Eq. 
(5.13) that for the anomalous diffusion 



Ar = (r q ) 1/q a 






(5.20) 



Expression (5.20) teaches us that in our model diffusion is anomalously fast over t , 
since a < 2 ; this diffusion is also called superdiffusion. At the same time it retains 
classical scaling over B ( that is, given by Eqs. (5.17) or (5.19)). We also note that 



the obtained t - dependence is the typical superdiffusion law within the framework of 
fractional kinetics, see Refs. [pi]] , p2f , |T8[] , p3| . 



Basing on the numerical solution to the Langevin equations (2.1) we estimated 
numerically the moments M r (t; q, a) by averaging over 100 realizations, each consisting 
from 50000 time steps. In Fig. 4 the g-th root of the g-th moment (characteristic 
displacement Ar) is shown versus t in a log-log scale for the three values of q, and 
for the Levy index equal 1.2, see from bottom to top: q = 0.12, which is less than 
a, q = 1.1 which is nearly a, and q = 2.0 (variance), which is greater than a. The 
dashed lines have the slope 1/a, which is, in fact, the theoretical value of the diffusion 
exponent at q < a, see Eq. (5.20). At q < a the numerical curve is well fitted by the 
dotted line. At q > a theoretical value of the moment is infinite, and in numerical 
simulation the moment strongly fluctuates. 

In our numerical simulation the scaling given by Eq. (5.20) was also checked in 
more detail. The results are presented in Fig. 5. At the left, in Figure a), we show the 
1/B - dependence of the characteristic displacement Ar. The Levy index a is equal 
1.2, and q is equal 0.12. The values of Ar obtained in numerical simulation are shown 
by black points, which are well fitted by straight dotted line. This fact confirms that 
Ar is inversely proportional to the magnetic field. At the right in Figure b) we show 
the exponent /i in the relation Ar ~ t M versus 1/a. The black points indicate the result 
of numerical simulation, whereas the dotted line shows the straight line fi = 1/a. We 
conclude that the right figure confirms the obtained theoretical dependence Ar ~ t 1 ' a . 



6 Conclusion 

In this paper we propose fractional Fokker-Planck equation (FFPE) for the description 
of the motion of a charged particle in constant magnetic field and stochastic electric 
fields. The latter is assumed to be a white Levy noise. We also assume that the particle 
is influenced by a linear friction force. Such formulation is a natural generalization of 
the classical problem for the Brownian particle [Q. It allows us to consider in a 



simplest way the peculiarities of the motion stipulated by non-Gaussian Levy statistics 
of a random electric field. 

The main results are as follows: 

1. The general solution to FFPE for a charged particle in constant magnetic field 
is obtained. In case of the absence of magnetic field this solution lead to the results 
obtained previously. However, the general solution also allows us to study the opposite 
case of a strong magnetic field in detail. 

2. The properties of stationary states are studied for two- and three-dimensional 
motions. The velocity relaxation is studied in detail, and non-Maxwellian stationary 
states are found, for which the velocity PDF is a Levy stable distribution. 

The energy PDFs are obtained, which have power law tails. This circumstance 
leads to divergence of the energy mean. In the real experiments as well as in numerical 
simulation the divergence manifests itself in large fluctuations of the mean value with 
time. 

3. The superdiffusion of a charged particle across magnetic field is studied within 
the framework of FFPE. The fractional moments of space displacement are estimated, 
and anomalous dependence of characteristic displacement Ar versus time, Ar ~ t 1 /", 
a < 2 , is found. The typical displacement is inversely proportional to the magnetic 



field. Therefore the diffusion described by our FFPE demonstrate anomalous behavior 
with time and remains classical with respect to the magnetic field dependence. 

4. We carry out numerical simulation based on the solution to the Langevin equa- 
tions for a charged particle in constant magnetic and random electric fields. We study 
numerically the process of relaxation and stationary energy states for different Levy 
indexes as well as super diffusion process. The results of numerical simulation are in 
qualitative agreement with analytical estimates. 

In summary, we believe that fractional kinetics will be a useful complementary tool 
for understanding and description of variety of non-equilibrium phenomena in space 
and thermonuclear plasmas. More elaborated fractional kinetic equations based on 
more sophisticated Langevin equations can be constructed, which, in particular, lead 
to a finite variance of the displacement of a particle and to anomalous B - dependence. 
On the other hand, a consistent development of the theory of essentially non-Gaussian 
plasma fluctuations is also of interest. 



Acknowledgments 

This work is financed by INTAS Call 2000, Project No. 847. One of us (M.S.) was 
supported by the KBN grant 2P03D01417. 



References 

[1] Bouchaud J. - P., Georges A., Anomalous Diffusion in Disordered Media: Sta- 
tistical Mechanisms. Models and Physical Applications, Phys. Reports 195, Nos. 
4&5, 127-294 (1990). 

[2] Metzler R., Klafter J., The Random Walk's Guide to Anomalous Diffusion: A 
Fractional Dynamics Approach. Phys. Reports 339, No. 1, 1-77 (2000). 

[3] Mandelbrot B. B. The Fractal Geometry of Nature. New York: Freeman, 1982. 

[4] Feder J. Fractals. New York: Plenum Press, 1988. 

[5] Mandelbrot B. B., Ness van, Fractional Brownian Motions, Fractional Noises and 
Applications. SIAM Review, 10, No.4, 422-437 (1968). 

[6] Levy P. Theorie de 1 'Addition des Variables. Paris: Gauthier - Villiers, 1937. 

[7] Skorokhod A. V. Random Processes with Independent Increments. Moscow: 
Nauka, 1986 [Engl, transl. Dordrecht: Kluwer, 1991]. 

[8] Maejima M., On a Class of Self-Similar Processes, Zeitschrift fur Wahrschein- 
lichkeitstheorie und verwandte Gebiete 62, 235-245 (1983). 

[9] Samorodnitsky G., Taqqu M. S. Stable Non-Gaussian Random Processes. New 
York: Chapman & Hall, 1994. 

[10] Chechkin A. V., Gonchar V. Yu., A Model for Persistent Levy Motion. Physica A 
277, 312-326 (2000). 



[11] B. V. Gnedenko B. V., Kolmogorov A. N. Limit Distributions for Sums of In- 
dependent Random Variables. Moskva: Izd-vo tekhniko-teor. lit-ry, 1949 [Engl, 
transl.: Addison Wesley Reading, MA, 1954]. 

[12] Mantegna R. N., Stanley H. E., Stochastic Processes with Ultraslow Convergence 
to a Gaussian: The Truncated Levy Flight, Phys. Rev. Lett. 73, No.22, 2946-2949 
(1994). 

[13] Chechkin A. V., Gonchar V. Yu., Self and spurious multi-affinity of ordinary Levy 
motion, and pseudo- Gaussian relations, Chaos, Solitons and Fractals 11, 2379- 
2390 (2000). 

[14] Korsunoglu B., Brownian Motion in a Magnetic Field, Annals Phys., 17, 259-268 
(1962). 

[15] West B.J., Seshadri V., Linear systems with Levy fluctuations, Physica A113, 
203-216 (1982). 

[16] Peseckis F. E., Statistical dynamics of stable processes, Phys. Rev. A36, No. 2, 
892-902 (1987). 

[17] Chandrasekhar S., Stochastic Problems in Physics and Astronomy, Rev. Mod. 
Phys. 15, 1 (1943). 

[18] Chechkin A. V., Gonchar V. Yu., Linear Relaxation Processes Governed by Frac- 
tional Symmetric Kinetic Equations, Journ. Eksper. Theor. Physics 118, No. 3(9), 
730-748 (2000). 

[19] Samko S. G., Kilbas A. A., Marichev O. I. Fractional Integrals and Derivatives, 
Theory and Applications. Minsk: Nauka i technika, 1987 [Engl, transl. Amster- 
dam: Gordon and Breach, 1993]. 

[20] Klimontovich Yu. L. Statistical Physics. Moscow: Nauka, 1982 [Engl, transl. New 
York: Harwood Academic Publishers, 1986]. 

[21] Jespersen S., Metzler R., Fogedby H. C, Levy Flights in External Force Fields: 
Langevin and Fractional Fokker-Planck Equations, and their Solutions, Phys. Rev. 
E59, 2736 (1999). 

[22] Zolotarev V. M., Uchaikin V. V., Saenko V. V., Superdiffusion and Stable Laws, 
Journ. Eksper. Theor. Physics 114, 1411-1425 (1999). 

[23] Yanovsky V. V., Chechkin A. V., Schertzer D., Tur A. V., Levy anomalous diffu- 
sion and fractional Fokker-Planck equation, Physica A282, 13-34 (2000). 



Figure Captions 

Fig. 1. Numerical solution to the Langevin equations (2.1). At the left: typical 
dependences of the energy E versus discrete time t for a) a = 2.0 (Gaussian noise term 
in the Langevin equations), and c) a — 1.1. At the right: (v x — v y ) plane for b) a = 2.0, 
and d) a = 1.1. 

Fig. 2. Stationary energy PDFs are shown in linear scale (at the top) and in log-log 
scale (in the bottom), for two-dimensional motion (in the left) and three-dimensional 
motion (in the right), respectively. The PDF with a = 1.1 is shown by solid curves, 
the PDF with a = 2 is shown by dotted lines. The points indicate the PDF obtained 
in numerical simulation. 

Fig. 3. The results of numerical simulation based on the solution to the Langevin 
equations. The q— th moments of the energy vs t for different orders of q and for 
a = 1.6. From bottom to top: q = 0.12 < a/2; q = 0.8 = a/2; q = 2.0 > a/2. In the 
bottom figure the dotted vertical line indicates the velocity relaxation time r v , whereas 
the dotted horizontal line shows the analytical value of the moment. 

Fig. 4. The results of numerical simulation based on the solution to the Langevin 
equations. The g-th root of the g-th moment as a function of t in a log-log scale for 
different orders of q and for a = 1.2. From bottom to top: q = 0.12 < a;q = 1.2 = 
a;q = 2.0 > a. The tangent of a slope of dashed lines equals to 1/a. 

Fig. 5. At the left: the g-th root of the g-th moment versus 1/B. The black points, 
which are fitted by straight dotted line, show the result of numerical simulation. At 
the right: the power /i in the relation (r q ) ' q oc t^ versus a. The black points, which 
are fitted by straight dotted line result from numerical simulation. 
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